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We study the influence of an external magnetic field on the collective electronic excitations in metallic 
nanoparticles. While the usual surface plasmon corresponding to the collective oscillation of the electrons with 
respect to the ionic background persists in the direction parallel to the magnetic field, the components in the 
perpendicular plane are affected by the field and give rise to two collective modes with field-dependent frequen- 
cies, the surface magnetoplasmons. We analyze the decay of these collective excitations by their coupling to 
particle-hole excitations and determine how their lifetimes are modified by the magnetic field. In particular, we 
show that the lifetime of the usual surface plasmon is not modified by the magnetic field, while the lifetime of 
the two surface magnetoplasmons present a weak magnetic-field dependence. Optical spectroscopy experiments 
are suggested in which signatures of the surface magnetoplasmons may be observed. 
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I. INTRODUCTION 



Collective excitations in confined many-body systems are 
of great fundamental interest. Such excitations decay due to 
their coupling to other internal degrees of freedom of the sys- 
tem. This allows the study of quantum dissipation and de- 
coherence and thereby the transition between quantum and 
classical physics. Particularly well studied is the case of sur- 
face plasmon excitations in metallic nanoparticles. 1 " 4 These 
collective dipolar vibrations of the electronic center of mass 
with respect to the ionic background dominate the optical ab- 
sorption spectrum. The corresponding resonance linewidth 
gives indirect access to the lifetime of the surface plasmon. 
Moreover, pump-probe experiments allow one to follow the 
time evolution of the electron dynamics with a resolution of 
a few femtoseconds, 5-7 and the surface plasmon excited by 
the probe laser field plays a prominent role in the interpreta- 
tion of the pump-probe measurements. 8 Recently, this tech- 
nique was combined with magneto-optical Kerr effect mea- 
surements to follow the magnetization dynamics in super- 
paramagnetic transition-metal nanoparticles and in particular 
the ultrafast demagnetization resulting from the pump laser 
excitation. 9 This demagnetization is also observed in ferro- 
magnetic thin films, 10 and its explanation is still a matter of 
debate. 11 - 12 

The saturation magnetization in ferromagnetic materials 
plays the role of an effective external magnetic field that cou- 
ples to the orbital degrees of freedom. 3 It is therefore relevant 
to study the influence of this effective field on the collective 
resonances in nanoparticles. In this work, we focus on the 
generic problem of the role played by an external magnetic 
field on the collective resonances in spherical nonmagnetic 
(i.e., alkaline or noble-metal) nanoparticles. While the surface 
plasmon excitation with the dipole parallel to the magnetic 
field is not modified, the two plasmon modes with dipoles per- 
pendicular to the field evolve in two magnetic -field-dependent 
collective modes when the magnetic field is switched on (see 
Fig. 1). We term these collective excitations "surface magne- 
toplasmons." They are split in energy by an amount of the or- 
der of the cyclotron frequency co c . Similar excitations exist in 
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FIG. 1. (Color online) Sketch of the center-of-mass motion with fre- 
quency a> M for the surface plasmon collective mode parallel to the 
magnetic field B, and for the two surface magnetoplasmon modes 
with frequencies a> + and <y_ [cf. Eq. (17)], where it rotates coun- 
terclockwise and clockwise in the plane perpendicular to B, respec- 
tively. 



quasi-two-dimensional semiconductor quantum dots. 14 ' 15 In 
the case of quantum dots, the energy scale of the resonance 
energies is not in the optical but in the infrared range, and the 
coupling of the collective excitation to the other degrees of 
freedom is very weak. 16 ' 17 In contrast, in metallic nanoparti- 
cles, this coupling is quite strong and leads to the decay of the 
collective modes. It is the aim of the present paper to analyze 
in detail the consequences of such a coupling on the lifetime 
of the surface magnetoplasmons. As we show in the sequel, 
the coupling yields a non-negligible linewidth of the corre- 
sponding resonances which might limit their observability in 
experiments. 

In this work, we analyze the decay of the surface mag- 
netoplasmons in metallic nanoparticles which is caused by 
their coupling to particle-hole excitations. This process, 
called Landau damping, 3 is the dominant decay channel for 
intermediate-size nanoparticles with a radius a in the range 
1 nm < a < lOnm. ls We show that the surface magneto- 
plasmon lifetimes (whose inverses yield the linewidths of the 
resonances) increase linearly with the size of the nanopar- 
ticle, as it is the case for the well-studied surface plasmon 
lifetime. 419-24 For experimentally available magnetic fields, 
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the surface magnetoplasmon linewidths are of the same order 
as the linewidth y of the usual surface plasmon. Since typ- 
ically cj c <K y, it is very difficult to individually resolve the 
collective magnetoplasmon resonances in a direct absorption 
experiment. However, we propose an indirect way to detect 
the predicted magnetic-field-induced modification of the ab- 
sorption spectrum by means of optical spectroscopy experi- 
ments. When the electric field of the exciting laser is lin- 
early polarized parallel to the external magnetic field, only 
the usual surface plasmon is triggered. In contrast, when the 
electric field is perpendicular to the external magnetic field, 
only the surface magnetoplasmons are excited (see Fig. 1). 
Using circularly polarized light, the two surface magnetoplas- 
mon modes can be individually addressed. The difference in 
the resulting absorption spectra should be clearly detectable 
with experimentally available magnetic fields. 

Furthermore, we argue that our results should be applicable, 
at least qualitatively, for ferro- or superparamagnetic nanopar- 
ticles, where the saturation magnetization is similar to an ef- 
fective external magnetic field. 13 It might therefore be neces- 
sary to include the effect of the magnetic field on the electronic 
center-of-mass motion when one attempts to model the effect 
of laser light on the magnetization of such nanop articles. 

The paper is organized as follows: In Sec. II, we present our 
model and the relevant collective coordinates to describe the 
surface magnetoplasmon dynamics. We calculate their life- 
times which constitute the main result of the present paper 
in Sec. III. We suggest in Sec. IV optical absorption spec- 
troscopy experiments to observe the surface magnetoplasmon 
collective modes before we conclude in Sec. V. Several tech- 
nical issues are explained in the appendices. 



is harmonic with the Mie frequency <d m = ^Ne 2 /m e a 3 in- 
side the particle and Coulomb-like outside. In Eq. (2), ®{x) 
denotes the Heaviside step function. 

Assuming that the magnetic field points in the z direction, 
B = Be z (see Fig. 1), we can express the Hamiltonian (1) in 
the symmetric gauge A = (— y, x, Q)B/2 as 



^ + U(n) + ^(i zj + 2s zj ) 



■s 



2|r, - r,-| 



(3) 



with (jl> c = eB/m e c the cyclotron frequency, and l z and s z the z- 
components of angular momentum and spin, respectively. No- 
tice that in writing Eq. (3), we omitted the diamagnetic term 
ZZi m e(^l( x ^ + y>)/% which is quadratic in the magnetic field 
and of the order of Nm e a>la 2 . As the Mie frequency typically 
lies in the optical range, and ha> c /B - 0.12meV/T, we have 
(jj c <K a>M, even for the highest presently achievable static 
magnetic fields of several tens of teslas. Thus, the diamag- 
netic term represents a small correction to the single-particle 
confinement (2), which is proportional to (w c /wm) 2 - In the 
sequel we work to first order in the small parameter (l> c /u>m- 



A. Separation into collective and relative coordinates 

A monochromatic electric field of (complex) amplitude Eo 
and frequency a> can be used to excite the electronic system. 
In the long-wavelength limit, the coupling to such an electric 
field is described by the Hamiltonian 



f^rrEoe 1 ' 



(4) 



II. MICROSCOPIC MODELING OF SURFACE 
MAGNETOPLASMONS 

We consider a spherical nanoparticle in vacuum containing 
N valence electrons of charge —e and mass m e . The nanopar- 
ticle of radius a = r s N 1 ^ 3 with r s = (3/47rn e )^ 3 the Wigner- 
Seitz radius (n e is the electronic density) is subject to a ho- 
mogeneous magnetic field B = V x A, where A is the vector 
potential. Within the jellium approximation, which assumes 
a homogeneous positive background of charge +Ne (Ref. 2), 
the electronic Hamiltonian reads 25 



Nr. , 

i=\ L e ' 

N e 1 



(1) 



with g e =i 2 the electronic g-factor, /l/b = efi/2m e c the Bohr 
magneton, c the speed of light, and s the spin operator. The 
single-particle confinement 

U(r) = ^{f~ 3« 2 ) ®(a - r) - ^0(r - a) (2) 



This field only couples to the electronic center of mass. 
Hence, it is appealing to decompose the Hamiltonian (3) 
by introducing the electronic center-of-mass coordinate R = 
Y,i *i/N and its conjugate momentum P = £j p,-. 24 ' 26 The rela- 
tive coordinates are denoted by = r, - R and pj = p, - P/N. 
Introducing this new set of coordinates in Eq. (3) and assum- 
ing that the displacement of the center of mass is much smaller 
than the nanoparticle size (|R| <K a), we obtain the decompo- 
sition 



H - H cm + H rcl + H c + R 



(5) 



where H cm and H re \ are the center-of-mass and relative- 
coordinate Hamiltonians, respectively, and H c is the coupling 
between them. The Zeeman term 



H z = w c ^ s Zti 



(6) 



accounts for the spin degrees of freedom. The decomposition 
(5) is reminiscent of the well-studied case where the degree 
of freedom of interest (in our case, the electronic center of 
mass) is coupled to a large reservoir or environment (the rela- 
tive coordinates). 27 The interaction with the reservoir leads to 
the dissipation of the collective coordinate energy. The num- 
ber of degrees of freedom of the reservoir is proportional to 
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the number of electrons in the nanoparticle, and thus, it is not 
very large. However, it has been shown 28 29 that such a reser- 
voir is sufficient to constitute a well-defined environment for 
the collective excitation, provided the nanoparticle is not ex- 
tremely small. 

The center-of-mass Hamiltonian appearing in Eq. (5) can 
be decomposed in two parts as 



H" +H ± 

cm cm' 



(7) 



with the part depending on the collective coordinate Z parallel 
to the magnetic field 



_ P 2 Z Mo? 
cm ~ 2M 2 ' 



(8) 



and the transverse component 



p2 , p2 



2M 
Nm e and L z 



Ma> 2 M 



(x 2 + Y 2 ) 



Y 



(9) 



XP Y - YP X - The frequency w M 



Here M 

<j>m N oa t/N is slightly redshifted as compared to the bare 
Mie frequency by the spill-out effect, 2-4,24 where N out is the 
number of electrons outside the nanoparticle. In what follows, 
we approximate w M by oj m for simplicity. The Hamiltonian 
for the relative coordinates reads 



rel 



z 



z 



211 



and the coupling Hamiltonian 



Ho = J]R-[VU(^]1 =q , 



with R • VC7(r) 



,R • r/(r) and 



f(r) = ®(a-r) + (^j 0(r- 



a). 



(10) 



(11) 



(12) 



The parallel part (8) of the center-of-mass Hamiltonian 
describes a harmonic oscillator which is independent of the 
magnetic field, and whose classical motion is sketched in 
Fig. 1. Its lowest excitation corresponds to the usual sur- 
face plasmon. 23 ' 24 In contrast, the transverse part (9) includes 
the effect of the magnetic field. The two associated collec- 
tive modes are called surface magnetoplasmons. Their cor- 
responding classical orbits are sketched in Fig. 1. Thus, an 
electric field which is polarized parallel to the magnetic field 
will excite only the usual surface plasmon in the z direction. 
In contrast, when Eq ± B, only the surface magnetoplasmons 
will be excited. In the general case, when Eo has components 
perpendicular and parallel to B, both the surface plasmon and 
the surface magnetoplasmons will be excited. Therefore, the 
relative orientation of the laser polarization and the external 
magnetic field allow one to selectively excite the different col- 
lective modes. 

The Hamiltonian (9) is similar to the one encountered in 
the context of quasi-two-dimensional semiconductor quan- 
tum dots. 1415 The main difference is that, unlike Eq. (2), the 



single-particle confinement in quantum dots is well approxi- 
mated by a harmonic potential for all relevant r. Hence, due 
to Kohn's theorem, 14 ' 15 ' 30 center-of-mass and relative coordi- 
nates decouple. In metallic nanoparticles, the Coulomb tail of 
the single-particle confinement (2) leads to a non-negligible 
coupling Hamiltonian (11) and therefore to the decay of the 
surface (magneto)plasmon excitations. Another difference is 
that, in quasi-two-dimensional semiconductor quantum dots, 
the relatively small effective electronic mass renders the cy- 
clotron frequency of the order of the confining one, such that 
the diamagnetic term in the Hamiltonian (1) cannot be omitted 
in the context of quantum dots. 



B. Effect of the magnetic field on the center-of-mass oscillation 



The Hamiltonian (8) can be written as 



1 



HL = hu M Wb+ - 



in terms of the bosonic operators 



1 / Z . Pz^osc 

b = — — - — + i 



a/2 \ Cc ft I' 



V2~Uosc ft /' 



(13) 

(14a) 
(14b) 



with ^ osc = y/h/MajM being the oscillator length. 

The Hamiltonian H£ m of Eq. (9) can be diagonalized by 
means of Fock-Darwin states. 14 Introducing the new variable 
£ = (X + \Y)I V2 and its complex conjugate and the bosonic 
operators 



(15a) 
(15b) 
(15c) 
(15d) 



b + = 


vl 


( ^osc 




bl = 




'A. 

( ^osc 


-I -\ 

cosc d?) 


b_ = 




'A. 








k ^osc 


bi = 




t?_ 

k ^osc 


-/ -) 

osc ^r 



we can write 



tfcm = ^ [b\b + + ^ + ruo- [blb_ + ^ . (16) 

The frequencies of the two magnetoplasmon excitations 
read 14 



(17) 



The surface magnetoplasmon with the larger (smaller) fre- 
quency a>+ (u>-) rotates counterclockwise (clockwise) in the 
plane perpendicular to the magnetic field (see Fig. 1). No- 
tice that one has lo+ > w_ since the Lorentz force -NeR x B 
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increases (decreases) the strength of the confinement seen by 
the collective mode with frequency a> + (ai_)- 

As is well known, 1415 the effect of the magnetic field is thus 
to split the usual surface plasmon in the plane perpendicular 
to the field axis in two surface magnetoplasmon excitations 
whose frequencies are separated by the cyclotron frequency 
u> c . However, the experimental observation of the two surface 
magnetoplasmon modes will be limited by the linewidths of 
these collective excitations. We address this important issue 
in Sec. III. 



C. Mean-field approximation for the environment 

The Hamiltonian for the relative coordinates (10) contains 
the electron-electron interactions. Assuming that the full cor- 
relations are not crucial for the present problem, we treat the 
interactions on a mean-field level. Numerical calculations us- 
ing the local density approximation and performed in the ab- 
sence of a magnetic field 23 suggest that the self-consistent po- 
tential treating the interactions on a mean-field level can be 
approximated by 



V(r) = V ®(r - a), 



(18) 



a spherical square well of depth Vo = sf + W, where ep and 
W are the Fermi energy and the work function of the consid- 
ered nanoparticle, respectively. We expect that the square well 
shape of the effective potential remains a good approximation 
for w c <s com- 31 We thus write, at mean-field level, 



#rel - 2 



2m e ' 2 z - 



(19) 



In what follows, we rewrite the above Hamiltonian together 
with the Zeeman Hamiltonian (6) in second-quantized form, 32 



#rel + H Z = 2^ Eao-C ] aa -C a 



(20) 



where the operators c acr {c aa ) create (annihilate) one-body 
eigenstates \a,&) with eigenenergies s acr = s a - crS<y c /2 in 
the mean-field potential V(r). Here a is a shorthand nota- 
tion for the orbital quantum numbers (n a , l a , m„), with n a , l a , 
and m a being the principal, azimuthal, and magnetic quantum 
numbers, respectively, while cr — +1 (-1) corresponds to spin 
up (down). It is important to realize that the Hamiltonians 
of Eqs. (19) and (20) do contain the electron-electron interac- 
tions (at mean-field level) and thus that our subsequent results 
concerning the surface magnetoplasmon lifetimes incorporate 
the crucial role played by the electronic interactions. 



D. Coupling of the center of mass to the environment 

Under our mean-field assumption presented in Sec. II C and 
using Eqs. (14) and (15), the coupling Hamiltonian (11) can 



be written as 

with the matrix elements 



<%e = (a,cr\zf(r)\/3,cr), 



(22a) 



€p ± = <r|(x ± iy)/(r)|0, cr), (22b) 

V2 

which contain the function f(r) of Eq. (12). The spherical 
symmetry of the wave functions associated with the Hamil- 
tonian (20) allows us to decompose the matrix elements (22) 
into angular and radial parts, 



dap ~ ^lj,h ^( e crer, s /3<t), 



(23a) 
(23b) 



The angular parts can be expressed in terms of Wigner 3 j 
symbols 33 as 



I la lp 1\| la h 1 



0) \—m a mp 



la lp 1|| la lp 1 



0/\— m a mp ±1/' 



(24a) 



(24b) 



with the selection rules l a = lp ± 1 for Eqs. (24), and m a = mp 
and m a = mp + 1 for Eqs. (24a) and (24b), respectively. The 
radial matrix element can be approximated by 21 



^(*W, £p<r) = 



cr 

m e a (Sacr ~ Spa-} 



2 ' 



(25) 



an expression that is obtained under the assumption of infinite 
confinement Vo — > oo in Eq. (18). 34 



III. SURFACE MAGNETOPLASMON DECAY RATES 



A. Landau damping of the surface magnetoplasmon collective 

modes 



For nanoparticles having intermediate sizes of the order 
of a few nanometers, the Landau damping 3 is the dominat- 
ing decay channel for the collective plasmon modes. 18 The 
surface (magneto )plasmons decay by producing particle-hole 
pairs in the electronic environment whose energies correspond 
to the ones of the collective excitations. Treating the cou- 
pling Hamiltonian (21) as a perturbation, the corresponding 
surface plasmon and surface magnetoplasmon decay rates, 
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whose inverses yield the lifetimes of these collective excita- 
tions, can be obtained from Fermi's golden rule. Assuming 
zero temperature, 35 we have the decay rates 23 ' 24 

T = j (mJosc^u) 2 \dph\ 2 $ { naj M ~ s p<t + s ha) (26) 



plUT 



and 



7± = ^ (we^osc^) 2 ^ \dp h J 2 6 {%0) ± - S pa - + E h(T ) (27) 



phcr 



for the surface plasmon and surface magnetoplasmons, re- 
spectively. Here, p and h denote particle and hole states with 
energies e pcr > sp and e/ lcr < sp. 

Using the expressions (23) for the coupling matrix ele- 
ments and the appropriate selection rules that are contained 
in Eq. (24), we obtain to first order in w c /wm <k 1 



y± = 



27T / 2/ ? osc 

fi \ a 

2.7T I 2,d osc 



(28a) 
(28b) 



where 



= f des(e -^ M )^a„,(e) 

(s) 



(29a) 



r 



-£p+tiUI ± 

max (sp,ti(i) ± ) 



X 



(^y+S ) Sl+l,rrer 1 ~ ^± ) 



(29b) 



Here Qi jn {s) is the density of states for fixed angular mo- 
mentum / and magnetic quantum number m. 22 For metals the 
Fermi wave vector kp ^ 10 run -1 , such that kpa » 1 for 
nanoparticles having a radius of more than a nanometer. Thus, 
we can resort to the semiclassical approximation 36 ' 37 to eval- 
uate the density of states. To leading order in ft, we have 23 ' 38 



, , y/2m e a 2 (s - huj c m/2)/H 2 - (I + 1/2) 2 

Ql,m (S) = 7T-, Z 7^7 . (30) 

which allows to evaluate Eq. (29) in the semiclassical limit 
(for details, see Appendix A), leading to 



r 



3 up 
47 



g(y) 



(31) 



for the surface plasmon decay rate. Here, Vp denotes the Fermi 
velocity, and v = Tium/ef- The function g{v) given in Ap- 
pendix B [cf. Eq. (B2)] is a monotonously decreasing func- 
tion. It is remarkable that the result of Eq. (31), which has 



been obtained to first order in the small parameter w c /wm, 
does not present any magnetic field dependence. We can thus 
conclude that for realistic field strengths the surface plasmon 
is not influenced by the external magnetic field. The result 
of Eq. (31) and its well-known 1 ja size dependence has been 
first obtained by Kawabata and Kubo 19 and refined later by 
many authors . 4 ' 20 ~ 24 

For the two surface magnetoplasmon collective modes, we 
find for cj c <K (wm,£f/S) (see Appendix A for details) the 
decay rates 



3 up 



cl> c 1 3 v , ^ 

(v)t — \^g(v)--g'(v)-h(v) i 
Li> M \2 2 j 



(32) 



which constitute the main results of our work. Here g'(v) de- 
notes the derivative of the function g(v) given in Eq. (B2), and 
the function h{v) is given in Appendix B [cf. Eq. (B4)]. The 
surface magnetoplasmon decay rates (32) present the same 
size dependence as that of the surface plasmon, Eq. (31), and 
obviously verify y±-y for zero magnetic field (u> c = 0). It is 
important to realize that the function 3g(v)/2 - vg'(v)/2 - h(v) 
entering the result of Eq. (32) is positive for v < 4.51, i.e., it is 
positive for any realistic values of the parameter v = Swm/sf. 
Thus, we conclude that y + (y_) decreases (increases) linearly 
for increasing co c (i.e., increasing magnetic field). This is 
due to the fact that the frequency co+ (<y_) of the correspond- 
ing mode shows the opposite behavior, i.e., it increases (de- 
creases) with u> c [cf. Eq. (17)]. The behavior of the decay 
rates y ± of Eq. (32) as a function of the magnetic field can 
easily be explained as follows: The typical dipole matrix ele- 
ments between particle and hole states separated by an energy 
As p h — e p - sy, entering the Fermi golden rule (27) scale as 
dph± ~ ^l^- e ph> w hile the density of particle-hole states ful- 
filling the selection rules dictated by the angular part (24b) 
of the matrix elements d^ h± is linear in Ae p /,. Since energy 
conservation requires that As p h = ti(±> ± , this argument, valid 
in the limit hu> c <K fiaiM <k sp (for details, see Ref. 28), 
where g(Q) = 1, g'(0) - 0, and h(Q) = 0, yields the scaling 
y ± ~ l/co 3 ± ~ 1 + 3w c /2w M , consistent with Eq. (32). 



B. Extension to dielectric environments and noble-metal and 
ferromagnetic nanoparticles 

In the case of alkali clusters in vacuum, the result of 
Eq. (31) for the surface plasmon decay rate agrees quantita- 
tively with numerical calculations using the time-dependent 
local density approximation as well as with experiments. We 
take this as a strong indication that the surface magnetoplas- 
mon decay rate (32) is quantitatively valid for alkaline clusters 
in vacuum as well. When one considers alkaline nanopar- 
ticles in a dielectric environment or noble-metal nanoparti- 
cles (in vacuum or in an embedding matrix), one has to take 
into account the steepness of the effective mean-field poten- 
tial V(r) entering Eq. (19) in the derivation of Eq. (31) to 
get quantitative agreement with numerical calculations as well 
as with experiments. 23 The resulting linewidth of the surface 
plasmon excitation decreases as the dielectric constants of the 
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metal and of the dielectric environment increase. This is also 
the case for the linewidths of the surface magnetoplasmons. 
Hence, their observation will be facilitated by a dielectric en- 
vironment. 

Furthermore, our results for the decay rates of the collective 
electronic excitations should, at least qualitatively, be applica- 
ble to ferromagnetic nanoparticles, provided one adds to the 
external magnetic field B the internal magnetic field 47rM s , 13 
with M s being the saturation magnetization. This internal field 
couples to the orbital degrees of freedom and is of the order of 
2T for bulk transition ferromagnets. Thus, in saturated ferro- 
magnetic nanoparticles, the surface magnetoplasmons might 
exist even in the absence of an external magnetic field. This 
possible extension of our considerations toward ferromagnetic 
nanoparticles might be important for the analysis of the mag- 
netization dynamics in such systems 912 since the surface mag- 
netoplasmons create electromagnetic fields inside the particle 
that might affect the magnetic moments responsible for the 
magnetism in these nanoparticles. 



IV. EXPERIMENTAL DETECTION OF SURFACE 
MAGNETOPLASMONS 

A. Absorption profiles 

For the nanoparticle radii of a few nanometer that we 
are considering, the extinction spectrum is dominated by 
absorption. 439 Assuming a Lorentzian profile for the line 
shapes of the collective modes, the (normalized) absorption 
cross section of photons with frequency a> is given by 



CT||(w) = 



y/2n 



(to - u M f + (y/2) 2 



(33) 



when the illuminating electric field is polarized parallel to the 
magnetic field [Eo = EqC z [| B, cf. Eq. (4)]. In that case, 
only the usual magnetic-field independent surface plasmon is 
excited. 

In the case where Eo _L B, the two collective surface mag- 
netoplasmon modes are triggered. Using linearly polarized 
light with Eq = £o(cos ip e A + sin if e v ), both magnetoplasmons 
are excited, and the absorption cross section reads 



o- ± (co) = - [cr + (a>) + cr-(co)] . 



with 



and 



o-+(co) = 



cr-(u>) = 



y+/2n 



(w- w+ ) 2 + ( r+ /2) 2 
y-/2n 

(w _ w _)2 + (r _ /2) 2- 



(34) 



(35) 



(36) 



The absorption profile <r ± , scaled with the maximum 2/ny of 
0"[|, is shown in Fig. 2 as a function of the photon frequency 
a> and for increasing cyclotron frequency a> c for the case of a 
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FIG. 2. (Color online) Absorption cross section cr ± for linearly po- 
larized light with electric field perpendicular to the external mag- 
netic field [Eq. (34)] as a function of the photon frequency a> and 
the cyclotron frequency ai c , both scaled with the Mie frequency 
w M . The chosen parameters correspond to a sodium nanoparticle 
(r s = 3.93«o, where ao is the Bohr radius) of radius a = 10 nm, 
having Tia> M = 3.5 eV. The horizontal white lines correspond to 
magnetic fields B = 10 T (dash-dotted line), B = 20 T (dotted line), 
B = 30 T (dashed line), and B = 40 T (solid line). The black lines 
indicate the surface magnetoplasmon resonance frequencies, (left 
line) and a> + (right line) [see Eq. (17)]. 



sodium nanoparticle of radius a = lOnm having a collective 
surface plasmon resonance at Sum = 3.5 eV. In the regime 
of currently available static magnetic fields, up to about 40 T 
(horizontal white solid line in Fig. 2), the two magnetoplas- 
mon modes (indicated by black lines in Fig. 2) are separated 
by a frequency a> c which is smaller than the surface magneto- 
plasmon linewidths y+ and y_ [Eq. (32)]. These linewidths are 
of the same order as that of the surface plasmon y [Eq. (31)]. 
Hence, the linewidths are larger than the separation of the res- 
onances such that the two modes may not be directly resolved 
in an absorption spectrum. 41 It should be noted that even us- 
ing pulsed magnetic fields that can reach up to 70 T would not 
allow one to clearly separate the two modes. As can be seen 
from Fig. 2, only for unrealistically large w c /wm > 1 % do the 
two collective excitations exhibit separate maxima. 

A way to individually address the surface magnetoplas- 
mons is to use circularly polarized light. In the case where 
E = E (e x + ie v )/ V2 [E = E (e x - ie y )/ V2], one can se- 
lectively excite the surface magnetoplasmon with frequency 
u> + [a>_], and the absorption spectrum is given by Eq. (35) 
[Eq. (36)]. But for currently experimentally achievable mag- 
netic fields, the displacement in frequency, as well as the mod- 
ification of the linewidth of these peaks, is so small that it 
might not be observable. 

In addition, for an ensemble of nanoparticles, the inho- 
mogeneous broadening due to their dispersion in size fur- 
ther increases the total linewidth. Thus, the direct and un- 
ambiguous observation of two distinct resolved peaks in the 
absorption spectrum of an ensemble of nanoparticles using 
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FIG. 3. (Color online) Differential absorption cross section Ao"ij n /cr 
of Eq. (37) as a function of the photon frequency to. The parameters 
are the same as in Fig. 2. The values of the magnetic field used in 
this figure correspond to the horizontal white lines in Fig. 2. 
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FIG. 4. (Color online) Differential absorption cross section Acr cir /<x 
of Eq. (38) as a function of the photon frequency to. The parameters 
are the same as in Figs. 2 and 3. 



linearly polarized light might not be possible. Moreover, us- 
ing circularly polarized light which enables one to excite a 
single surface magnetoplasmon mode may not present a suf- 
ficient magnetic -field-dependent behavior in the absorption 
cross section. However, we suggest in the following two dif- 
ferential measurements where one could experimentally iden- 
tify the effect of the magnetic field on the collective excita- 
tions and the presence of the surface magnetoplasmons, one 
using linearly polarized light (Sec. IV B) and the other, which 
is more efficient, using circularly polarized light (Sec. IV C). 
The effect of inhomogeneous broadening on these proposed 
differential measurements is discussed in Sec. IV D. 



B. Detection with linearly polarized light 

The relative differential absorption cross section 

given by the difference between the absorption cross sections 
for the electric field linearly polarized parallel [Eq. (33)] and 
perpendicular [Eq. (34)] to the external magnetic field is de- 
picted in Fig. 3 for the values of the magnetic field indicated 
by the white lines in Fig. 2 (the parameters are for the same 
sodium nanoparticles of radius a = 10 nm). It can be seen 
in Fig. 3 that the differential absorption A<xi in /cr increases for 
increasing magnetic field strength, yielding a clear signature 
of the existence of the surface magnetoplasmon excitations. 
The increase of the total linewidth of the absorption reso- 
nance caused by the magnetic field via the modified width and 
the splitting of the surface magnetoplasmons leads to positive 
values of AtT\i n /cr in the center of the resonance and to neg- 
ative values in the tails. In magnetic fields of less than 10 T, 
Ao"ii n /cr is very small, at least for the case of sodium nanopar- 
ticles. However, in a field of 40 T (solid lines in Figs. 2 and 
3), of the order of the highest static field available in present- 



day high magnetic field laboratories, the differential absorp- 
tion Acr lin /cr becomes noticeable. Furthermore, magnetic field 
pulses of 70 T are available with durations of the order of mil- 
liseconds that should be sufficient to measure the absorption 
cross section. For the latter value, we expect a relative differ- 
ential absorption of about 6 % that should be detectable (not 
shown in Fig. 3). 



C. Detection with circularly polarized light 

The relative differential absorption cross section 

= (38 ) 
O" <x + + (T- 

given by the difference between absorption cross sections for 
the two circular polarizations of the electric field, Eqs. (35) 
and (36), is shown in Fig. 4 for the same parameters as in 
Figs. 2 and 3. It is positive (negative) for a> > w M (a> < oj m ) 
due to the fact that the resonance frequency co + (<w_) of the 
surface magnetoplasmon mode "+" ("-") is larger (smaller) 
than the surface plasmon frequency ojm [cf. Eq. (17)]. As 
one can see from Fig. 4, the differential absorption cross sec- 
tion Acr c j r /Acr of Eq. (38) is much more pronounced than 
the one of Eq. (37) using linearly polarized light (compare 
with Fig. 3). Moreover, the frequency range where Acr^/cr 
presents noticeable values is much broader. In the case of 
circularly polarized light, Acr cir /cr is already of the order of 
several percent at a magnetic field of 10 T (green dash-dotted 
line in Fig. 4) and can reach up to about 15 % in a static field 
of 40 T (black solid line). For a pulsed magnetic field of 70 T, 
we obtain a maximal relative differential cross section of about 
25 % (not shown in the figure). It thus seems even more effi- 
cient to use circularly polarized light rather than linearly po- 
larized light to detect the surface magnetoplasmon modes in 
metallic nanoparticles. 
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FIG. 5. (Color online) Differential absorption cross section A<x cir /cr 
of Eq. (38) as a function of the photon frequency a> in a magnetic 
field B = 10 T for an ensemble of sodium nanoparticles with mean 
radius a = 10 nm and for increasing values of the inhomogeneous 
linewidth y m . The inset is the same as the main figure for linearly 
polarized light, Eq. (37). 



D. Effect of inhomogeneous broadening 

In Sees. IV A, IV B, and IV C, we have implicitly assumed 
that the nanoparticles irradiated by the laser light do not 
present any size dispersion. Although this is a valid approx- 
imation for experiments on single clusters, 42 ^ 7 most experi- 
ments are done on ensembles of metallic clusters where the 
inhomogeneous broadening resulting from the size depen- 
dence of the resonance frequency masks the homogeneous 
linewidth. 48 " 50 

In this section, we address the role of inhomogeneous 
broadening on the experimental detection of the two surface 
magnetoplasmon excitations. To this end, in the absorption 
cross sections of Eqs. (33), (35), and (36), we phenomeno- 
logically add an inhomogeneous linewidth y m resulting from 
the size dispersion of the ensemble of nanoparticles to the in- 
trinsic linewidths of the surface (magneto )plasmon collective 
modes. The resulting differential absorption cross sections for 
sodium nanoparticles with mean radius a = lOnm in a mag- 
netic field of 10 T using circularly polarized light [cf. Eq. (38)] 
are presented in Fig. 5. Results for the same parameters using 
linearly polarized light [cf. Eq. (37)] are shown in the inset in 
Fig. 5. 

We quantify the inhomogeneous broadening y m by its rel- 
ative value with respect to the intrinsic linewidth of the col- 
lective modes at zero magnetic field, y [see Eq. (31)]. Given 
the size dependence of the resonance frequency of the sur- 
face plasmon excitation which can be estimated from numer- 
ical calculations within the time-dependent local density ap- 
proximation (see Fig. 7 in Ref. 24), and assuming a constant 
size distribution with dispersion Aa around the mean radius 
a, we find that Aa/a = 50% (Aa/a = 100%) corresponds to 
Tin/y = 0.5 (yin/y = 1) for the parameters of Fig. 5 (see blue 
dotted and red dashed lines). 

As one can see from Fig. 5, inhomogeneous broadening 



tends to weaken the differential absorption cross sections. In- 
deed, the maxima of the differential cross sections for cir- 
cularly (Fig. 5) and linearly polarized light (inset in Fig. 5) 
decrease for increasing inhomogeneous linewidth y m . For 
an inhomogeneous linewidth y m comparable to the intrinsic 
linewidth y, the maximal differential absorption cross section 
using circularly polarized light, Acr cir /cr, is still of about 2% 
in a magnetic field of 10 T (see red dashed line in Fig. 5). 
Such Acr c i r /cr should still be clearly measurable. 49 50 In the 
case of linearly polarized light, the maximal differential ab- 
sorption cross section is only about 0.03 % for y m = y (see 
red dashed line in the inset in Fig. 5). Such a small value in 
the difference between two extinction spectra on ensembles of 
nanoparticles seems to be difficult to measure. 49 ' 50 These re- 
sults point out that using circularly polarized light to obtain 
a clear-cut experimental detection of the two surface magne- 
toplasmon excitations in ensembles of metallic nanoparticles 
presenting a rather large size (and shape) dispersion is more 
appropriate than using linearly polarized light. 



V. CONCLUSION 

We have analyzed the role of an external magnetic field on 
collective excitations in metallic nanoparticles. As is the case 
in the context of semiconductor quantum dots, the magnetic 
field induces two new resonances in metallic nanoparticles, 
the surface magnetoplasmons, which can be excited when the 
polarization of the electric field has a component perpendicu- 
lar to the magnetic field. Our main result concerns the Landau 
damping linewidths of these collective modes that we have 
calculated. In particular, we have shown treating the electron- 
electron interactions within a mean-field approximation how 
the magnetic field modifies the absorption linewidths of the 
surface magnetoplasmon resonances. In all realistic cases, the 
linewidths are much larger than the splitting of the resonance 
energies such that one may not resolve them directly in a sin- 
gle absorption measurement. Nevertheless, if one changes the 
polarization of the electric field with respect to the magnetic 
field from parallel to perpendicular, a noticeable change in 
the absorption profiles should be detectable, at least in very 
strong magnetic fields, as they are currently available in high- 
field laboratories. Using circularly polarized light, which en- 
ables one to selectively excite the surface magnetoplasmon 
modes, leads to even larger values of the differential cross sec- 
tion when one changes the polarization from right-handed to 
left-handed. The proposed differential measurements are ex- 
pected to allow a detection of the surface magnetoplasmon 
modes even in the presence of an inhomogeneous broadening 
of the resonances that results from the size and shape disper- 
sion for ensembles of nanoparticles in matrices. 

Our results can be extended to ferromagnetic nanoparticles, 
where the internal magnetic field (i.e., the saturation magne- 
tization) couples to the orbital degrees of freedom. In addi- 
tion, when the ground-state magnetization in nanoparticles is 
nonzero, the collective electronic excitations couple to spin- 
dependent excitations in the nanoparticle. 51 This is the case 
in nanoparticles with an open electronic shell, particularly 
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for the case of ferromagnetic nanoparticles. It then becomes 
possible to affect the magnetization indirectly by exciting the 
charge degrees of freedom, and the influence of the magnetic 
field on the latter might become important. 
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Appendix A: Semi classical calculation of the surface plasmon 
and surface magnetoplasmon linewidths 

In this Appendix, we present the details of our semiclassical 
calculation of the surface (magneto)plasmon linewidths dis- 
cussed in Sec. Ill A. Introducing the notation so = h 2 /2m e a 2 
and k — fui> c /2so, we expand the semiclassical density of 
states for fixed / and m [Eq. (30)] for e » fia> c to obtain 52 



and 



£>/,m(e) - Qifiis) [1 + Km f i{e)} , 



where 



i 



(Al) 



(A2) 



s/e 2 e / £() -(/ +1/2) 2 ■ 

With Eq. (Al) and using the expression of the angular matrix 
elements (24), the functions F and F± [cf. Eq. (29)] that enter 
the surface plasmon and surface magnetoplasmon linewidths 
[cf. Eq. (28)] read 

F = - \ dss(s- ftu> M ) 2^ Qi,o(e) 

J max (sp,%0m) ; 

X [(I + l)£>/+i,o(e - Sw M ) + lQi-\,o(s - Sw M )] , (A3a) 



and 



F ± 



1 p*** \— i 

- dss(s-hu) ± ) 2^eifi(s) 

3 xj max (ep fito± ) / 

x |(Z + l)e/+i,o(e - fr(jL)±) 1 + Kfi+i(e - tioj ± ) 
1 + k/i-i(e - ho> ± ) 



+ lgi-\fl(e - tiu> ± ) 
l+l 



K\fi{s - fki) ± ) + fi-\{e - ho> ± )j 



(A3b) 



Consistently with the semiclassical approximation (high- 
energy limit), we now assume that / » 1 (i.e., / I ± 1) and 
approximate in Eq. (A3) the summation over / by an integral. 
With Eqs. (30) and (A2), we obtain 

F = de Al21 



(A4a) 



Jmax(£ F; ftw ± ) JO 



■ V(£-S«±)/£o 



dZ 2/ 



so 



X < 1 ± K 



So 

1 



So 



1 



2 (e - tu±> ± )/so - I 2 s-Tiu ± 



(A4b) 



Expanding the above expressions for o) c <K one finds with 
Eq. (28) our final results for the surface (magneto)plasmon 
linewidths, Eqs. (31) and (32). 



Appendix B: Auxiliary functions for the surface plasmon and 
surface magnetoplasmon linewidths 

The function g(v) entering the surface plasmon linewidth 
(31) is defined as 



g(v)=- f dx f dyy/(x-y)(x-y-v). (Bl) 

v Jmax(l,v) JO 

The double integral can easily be evaluated, and one finds 21 

g (v)= _L[ ( l +v) 3/ 2 _ (1 _ v) 3/2] 



3v 
v 

4 
v 



+ - ( Vl + v - Vl - v - vln v) 

(i + £)in(vnrp-i) 
-(i-5)m(i-vr^)" 

for v < 1 and 

*(v) = ^(1 + v) V2 + - ( VTT^ - In y) 



(B2a) 



+ 



(l + 01n(VTT^-l)-^ln Vv 



(B2b) 



for v > 1. The function g(v) is shown in Fig. 6 (red dashed 
line). Its asymptotic behaviors are g(v) 1 + v 2 [ln(v/4) - 
l/6]/4 and g(v) m 8/15 yfv for v <«; 1 and v » 1, respec- 
tively. Its derivative g'(v) entering the surface magnetoplas- 
mon linewidth (32) is shown as a blue dotted line. 

The function h(v) entering the surface magnetoplasmon 
linewidths (32) is defined as 



h( 



pl+y px—v 

v) = | dx I dy y/(x - y)(x - y - v) 

Jmax(l,v) JO 
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FIG. 6. (Color online) The function g(v) of Eq. (B2) (red dashed line), its derivative g'(v) (blue dotted line), and the function h(v) of Eq. (B4) 
(black solid line). 



Its explicit expression reads 
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(B4a) 



(B4b) 



for v > 1 . It involves the dilogarithmic function 

ln(l -f) 



Li 2 (z) 



df ■ 



I 

The asymptotic behaviors of the function /i(v) for v <K 1 and 
v » 1 read /z(v) ^ -v 2 [ln(v/4) + l]/4 and h(v) 2 yv/9, 
respectively. The function h(v) is plotted in Fig. 6 (black solid 
(B5) line). 
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